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Abstract 

We formulate the conformal bootstrap approach to four-fermion theory at 
its strong coupling fixed point in dimensions 2 < d < 4. We present a solution 
of the bootstrap equations in the five-vertex approximation. We show that the 
bootstrap approach gives a particularly simple way to obtain next to leading 
order corrections to critical exponents in the large N expansion and present the 
values of the anomalous dimensions of the fermion field ip and the composite 
iptp to order 1/N 2 . 
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Conformal symmetry describes the critical behavior of quantum field theories or 
statistical systems at points of second order phase transition. It is widely used in two 
dimensions where it provides exact solutions to a variety of non-trivial interacting 
field theories. In this Letter we shall demonstrate that it can also be useful in higher 
dimensions. If a field theory has a fixed point of the renormalization group flow, the 
dynamics at that point is conformally invariant. In any spacetime dimensions, confor- 
mal symmetry determines the form of the two and three-point Green functions of the 
field theory up to a few constants, the scaling dimension of the field operators and the 
value of the coupling constant at the fixed point. This information, when combined 
with Dyson-Schwinger equations can provide a powerful tool for the analysis of field 
theories. 

In fact, in a field theory with only three-point coupling constants, the Dyson- 
Schwinger equations which relate two and three point functions could in principle be 
used to determine all anomalous dimensions in terms of the coupling constants, which 
are then left to be determined by higher order equations. Unfortunately, this direct 
approach is plagued by technical difficulties involving indeterminate expressions. This 
problem was partially solved by Parisi |]J who showed how to use the information 
of the lower order Dyson-Schwinger equations to derive bootstrap equations for the 
effective coupling constants. The resulting equations can not be solved exactly, but 
can be used to find interesting approximate solutions of a field theory. 

In this Letter, we shall consider the example of a fermionic field theory with 
a four-fermion interaction, — 4p {j2iL\ ^i^ij m spacetime dimensions 2 < d < 4. 
By introducing an auxiliary scalar field <p = <^Si=i'0i'0i (with no kinetic term), 
this is equivalent to a model with a Yukawa vertex —X(j)Y,f=i'4'i' l Pi + 2 /2- The 
renormalization group flow in this model has an infrared stable fixed point at A = 
and there is a long-standing conjecture that it also has an ultraviolet stable fixed point 
at some A* > 0. This fixed point is seen both in the large N and in the e expansion. 
Also, as d — > 2 + , A* — > to obtain the asymptotic freedom of the Gross-Neveu model. 
Recently, renormalizability of the large N expansion of d = 3 four- Fermi theory has 
been proven using the rigorous renormalization group technique of constructive field 
theory [g]. 
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In this Letter we shall assume that A* exists and analyze the resulting conformal 
field theory. We find a solution of the bootstrap equations in the three and five 
vertex approximation. In particular, we show that these give an easy way to find the 
anomalous dimensions of ip and in the large N limit to order 1/N 2 . 

Under rescaling of the coordinate, x — > px, the fields transform as 

iP'(px)=p- l iP(x), ft(px)=p- l $(x), <j>'(px)=p- b <j>(x), (1) 

with I and b the scaling dimensions of the fermion (anti-fermions) and the boson, 
respectively. Under a special conformal transformation parameterized with a constant 
vector t^, 

rp __, - x v + t ^ x (n\ 
»~ l + 2t-x + t 2 x 2 ' [ ) 

where a x = cr t (x) = Ij^rl 1 ^ = 1 + 2t ■ x + t 2 x 2 , they transform as 

il>(x) -> i/>'(x') =a l x - 1/2 {l + tx)^{x) , (3) 
$(x) -> i)'{x') =a l - l ' 2 ^(x)(l+xi) , (4) 
- ^V) = (5) 

where x = ^^x^ (7^ are the Dirac Matrices). 

These, together with translation and Lorentz invariance, form the global finite 
dimensional conformal group. Conformal invariance determines the form of two-point 
correlation functions. In momentum space, 

g(p) = kp 2 y- h+i/2 , D^) = - 2 ( P 2 r h+i . (6) 

A convenient normalization has been chosen. (This can always be done using a finite 
rescaling of <j) and ip.) Moreover, conformal invariance fixes the three point vertex up 
to a dimensionless constant factor, the effective coupling A, 

v( \ — \ N ^ r^h k+f>i h+f>2 1 (7) 

l\Pl,P2) N(p)N{l — 6/2) J 7T h [(k + Pl ) 2 ] 6/2+1/2 [{k + p 2 ) 2 ] b /2+V2 [fc2]I-6/2 ' ^> 

where N(t) = T( ^^ ', and 

7 = Z + 6/2 - /i , (8) 
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which is defined as the index of the vertex. By dimensional analysis, it is easy to 
check that the (anomalous) dimension of the vertex T in momentum units is —2 / ~f. 

Our task is to determine the scaling dimensions I and b, which describe the critical 
behavior of the model, and the critical coupling constant A = A*. We shall use the 
bootstrap equations. These are derived from the Dyson-Schwinger equations for the 
fermion and boson self-energies and for the Yukawa vertex. We present the one for 
the vertex first. Graphically, the bootstrap equation of the Yukawa vertex is ||[IJ] 




Fig. 1 The bootstrap equation for the Yukawa vertex. The dark spots represent 
the conformal Yukawa vertex, the solid and dashed lines represent the conformal 
fermion and boson propagators, respectively. 



It is worth noting that the contribution of the bare vertex which would normally 
appear in the right-hand-side of this equation is absent. This is a result of the infinite 
renormalization (and therefore zero of the multiplicative renormalization constant) 
which is necessary to make the theory finite. Moreover, it is required by the conformal 
ansatz, since the bare vertex does not have the conformally invariant form given in 
([7]) P|. The Dyson-Schwinger equation can be rewritten in a neat form by factoring 
the momentum dependence from both sides. To see this, let the momentum carried 
by the external boson field be zero. Then the l.h.s. is readily calculated. Performing 
the integration in (^) with pi = p 2 = p, we obtain 

nv,v) = ^- (9) 

The r.h.s. of the bootstrap vertex equation is an infinite series. By dimensional 
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analysis, each term is proportional to j^y. More precisely, the n-th term is 

A (2 " +1) /(2n + l)(M,&)^, (10) 

where the function /(2n+i)(^ I, b) is given by the diagrams with 2n + 1 conformal 
vertices. The arguments /, I, and b refer to the scaling dimensions of the external 
fermion, anti-fermion, and boson, respectively. The bootstrap equation for the vertex 
is 

1 = \ 2 J(l,l,b;\*), (11) 

oo 

f(l, I, b; K) = A^ n-1 )/( 2n+ i)(Z, I, b) . (12) 

n=l 

We call f(l, I, b; A*) the vertex function. 

There are two similar Dyson-Schwinger equations for the fermion and boson self- 
energies, respectively, 



S = — »- (a) 




Fig. 2 Dyson-Schwinger equations for the fermion (a) and boson (b) self-energies. 
The vertex without a dark spot stands for the bare Yukawa vertex, associated 
with the bare coupling Ao- 

As one might expect for primitively divergent vertices in a scale invariant theory, 
these are ill-defined. The l.h.s. is finite. The r.h.s. is a diverging integral times a 
multiplicative renormalization constant and is the indeterminate • oo. A way to 
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resolve this ambiguity was suggested by Parisi JI]]. It involves taking derivatives of 
each side by external momenta and by the scale dimensions of the operators. The 
result is the bootstrap equations for the vertex function. In the present model, we 
obtain |§ 

Xi N{ 7 )N(d - l)N(l)N\b/2) df{V, I, b; A* 



1 



(4ir) h V(h) N 2 (b)N(l -b/2) dl'/2 

Xi N( 1 )N 2 (b/2)N(d-b) df{l,l,V;\, 



i'=i 



iVTrl 



(13) 
(14) 



b'=b 



(47r) h T(h) N(b)N(l - b/2) dV/2 

where N(r) is defined after ([7|) and N(t) = • Also, N is the number of 

fermion species and Trl is the dimension of the Dirac Matrices. 

It is remarkable that, in all three equations (|T3|), (|14|) and (11]), the only unknown 



is the vertex function f(l,l, 6; A*), which contains all the necessary information to 
determine the three parameters. In particular, there is no need to calculate any 
Feynman diagrams for the self-energies. 

Unfortunately, calculating the vertex function f(l,l, b; A*) is a difficult problem, 
as it involves an infinte number of Feynman diagrams. On the other hand ([12]) is a 
power series in the coupling constant and if the critical coupling is small, the first 
few terms are expected to give good approximate results. One limit in which A* is 
small is the limit of large number of fermion species, N. There, A* ~ 1/N + ■ • ■ and, 
as we shall demonstrate below, further approximation to the bootstrap equations can 
be solved analytically to order 1/N 2 . Quantities which are also small in the large iV 
limit are the anomalous dimensions defined by 

7^ = l-h + 1/2, (15) 

7<a = 6-1; (16) 

which are both of order 1/N. In terms of them, the index of the vertex is 

1 = l + b/2-h = llp + l4> /2 . (17) 

The integral on the r.h.s. of the vertex (|7|) is related to the Appel function F 4 ||. 
It takes a simple form, when p\ < ^ip\= p 2 . In this case, 

A N(b/2)N(l) , . 

V (pi= ,, n=p) = — w ^lAl- ) . (18) 



It is easy to check that, at the lowest order in 7, ( pjj|) coincides with ((^). 

To solve the set of conformal bootstrap equations, we start with calculating the 
three-vertex correction to the vertex function. Using the conformal propagators and 
vertex, and setting the momentum carried by the external boson zero (i.e. p% = P2), 
the three-vertex diagram in Fig. 1 is 

d d k 



-A 



(27r) d 
A 3 1 



T(p,k) 



k 



;T(k,k)- 



k 



N(b/2)N(l) , 2 N(j + b) 



T72-r(fc,p) 



l(k-p) 2 ] h - 



(19) 



(p*)i (An) hV N(b)N(l-b/2y N(b)N(j) ' 
Above, we have used ( |18D for T(p, k), as the integration over the region of large 



internal momentum dominates. Then we have 

1 



fs(l,l,b) 
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+ ... 



(20) 



(47r) A r(/i)7 L ' h-1 

Substituting the first term of (^) into the set of conformal bootstrap equations 
TIP, (|13|), and (14), we obtain the anomalous dimenions of the fermion field ip and 
composite operator = ipip at the leading order 1/N 

1 F(2h - 1) sin(7r/i) 



7? 



,(!) 



ATTrl 7rr(/i - l)T(/i + 1) ' 
2(2/i-l) (1) 



"7^ 



(21) 
(22) 



which reproduce the known results, with 

2 (1) 1 (47r)*r(2/i - 1) sin^/i) 

7 (1) 



A^Trl 7iT(h) 
1 r(2/i - 1) sin(7r/V) 



(23) 
(24) 



iVTrl 7rr 2 (/t) 

Moreover, to the next to leading order, as we shall see below, the vertex function 



/ takes a same form as that in (|20|) . This implies to the order 



df(l',l,b; K) 



dl'/2 



df(l,l,b';\* 



i'=i 



db'/2 



(25) 



b'=b 



Then, with no need of further details of the vertex function f(l,l,b), by using the 



boostrap equations (0) and fll4j) , one can determine the fermion anomalous dimenion 



to order 1/N 2 . Dividing one of them by the other, we have 

1 N(l)N(d-l) 



(26) 



NTrl N(b)N(d - b) ' 
Expand the r.h.s. of the above equation over 7^ and 7^ (both ~ 1/N+...), we obtain 

1 

7^ = 7i 1} [l + \~ +^(2h-l)-m + rrctg(hrr)) + ...], (27) 

where -?/>(x) = r'(x)/r(x), and 7^, and 'j^ are given in (pi]) and fl2"2"|). 

To obtain the anomalous dimension of = to order 1/N 2 , one needs the vertex 
function / to the same order. Besides the second term in the three-vertex correction 
fl2"0|), it involves the leading order of the five- vertex correction, and also the leading 
order of the seven-vertex diagrams with a fermion loop. 

Notice that only one five- vertex diagram, depicted in Fig. 1, contributes to f 5 (l, I, b). 
(Another diagram which also contains five conformal vertices has a subdiagram with 
a fermion loop attached to three boson legs and is zero by parity symmetry). The 
Feynman integral (setting p\ = P2) is 

G(k + q- p)T(k + q-p, q)G(q)T(q,p)D(k - p)D(k - q) . (28) 
Performing the integral, to the leading logarithmic term, we obtain 

(29) 

To order 1/N 2 , the seven-vertex diagrams with a fermion loop attached to four 
boson lines contribute, as a fermion loop carries a factor N. This sort of seven-vertex 
diagrams are given in Fig. 3. 






A^A r*\ A^A 



Fig. 3 The diagrams of seven-vertex correction with a fermion loop. Each fermion 
loop carries a factor N. 



For f-j in the second order, to the leading logarithmic term again, we have 
NTtlT(2-h) 1 3r 3 (fe/3)r(h)(2fe/3-l) | ^ 



(4tt) 3/i (/i- l) 2 r(2/i-2) 7 L 8/ir 3 (2/i/3) 

jvrtir(2-/i) i | i r 3 (fe/3)r(fe) | 



(4tt) 3,1 (/i- l) 2 r(2/i-2) 7 L 2r 3 (2/i/3) 
for the first and the next two diagrams in Fig. 3, respectively. The vertex function is 

/(/, I, b, A) = f 3 (l, I, b) + X 2 J 5 (l, I, b) + Xif r (l, /,&) + ..., (32) 



where f 3 , f 5 , and f 7 , to the second order, are given in (|2Cf), (|29|) , and (|30| ) and fl31|) . 
This is a systematic expansion procedure. To calculate the next order, 1/iV 3 , for 
instance, requires the first three terms in the expansion of f 3 , the first two terms of 
/s, and the first two terms of (one with a factor N carried by a fermion loop and 
the other without), and the first term of fg with an explicit N factor. Similarly, as 
the multiplicative factors of the partial derivatives of /(/, I, b; A*) of equations flT3| ) 
and flTj] ) are functions of I and b, we expand these to the second order in 7^ and 7^ 
(or 7) as well. Substituting these expansions in the equation ( |i"3"D and (]TT|) , we obtain 

A* = -^(Tcctgi^ + ^h-l)- 4,(1)) + ...], (33) 

7 ^ - ^ + _ 1} - if^h) + ^(2* - 1) - V(l)) 

(A 2(1) ) 2 iVTri r(2 - h)r(/i) ,2/1 - 1 r 3 (V3)r(/>) q + 1 _ jj_ + ^ _ 



(47r) 2/i r(2/i-2)(/i- l) 2 v /i-l 2r 3 (2/i/3) v 2 V 2/i y 

(34) 

Then finally, 

70 = 2(7 - 7v) > (35) 
where 7 and 7^ are given in (|3~41 ) and (P7|). 

In particular, when d = 2h = 3 and Trl = 2, we have 

2 r 148 

7 * = 3^iV [1 + 9^iV + --- ] ' (36) 

16 . 76 + 27tt 2 . , , 

^ = -3^n [1 + ^6^t + - ] - (37) 
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The i/j anomalous dimension (7^,) has been computed to this order by conventional 
diagrammatic methods in ||. That calculation involves much more labor than the 
present one. Also, the present calculation gives the anomalous dimension of the 
composite operator = ipip which has not been computed elsewhere. Note that the 
result for 7^, is not the same as in ||. (The numerator of the second term quoted 
in P has 122 instead of 148.) We presently do not understand the source of this 
discrepancy. 

Here, we see the power of the bootstrap method in performing approximate cal- 
culations, such as the large N expansion to higher orders. It would be interesting 
to apply this approach to other field theories, such as gauge theory or theories with 
vector-like vertices. It would also be interesting to search for other, non-perturbative 
solutions in the present four-fermion, or in A0 4 theory [^]. 

Note added: After this work was submitted for publication, J A Gracey informed 
the authors that the number 122 was a misprint in ||, his result would be 112, instead. 
In a recent publication |§j he also calculated the exponent 7^ to order 1/iV 2 . 
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